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Abstract 

We present some analytical solutions to the Einstein equations, de- 
scribing radiating collapsing spheres in the diffusion approximation. 
Solutions allow for modeling physical reasonable situations. The tem- 
perature is calculated for each solution, using a hyperbolic transport 
equation, which permits to exhibit the influence of relaxational effects 
on the dynamics of the system. 
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1 Introduction 



Our purpose in this work is to provide some analytical solutions to Einstein 

equations, describing collapsing dissipative spheres in the difussion approx- 
imation. Such solutions may serve as test-bed for numerical relativity, and 
for probing cosmic censorship and hoop conjecture, among other important 
issues, and represent a natural tool to bring out the influence of dissipation 
on the evolution of a collapsing object. 

Analytical solutions although generally are found, either for too simplis- 
tic equations of state and/or under additional heuristic assumptions whose 
justification is usually uncertain, are more suitable for a general discussion 
than purely numerical solutions, which usually hinder to catch qualitative, 
aspects of the process. 

Therefore it seems useful to consider analytical models which are rela- 
tively simple to analyze but still contain some of the essential features of a 
realistic situation. 

Our endeavour is further justified by the following two considerations: 

• It is already an established fact, that gravitational collapse is a highly 
dissipative process (see [1], [2], [3] and references therein). This dissipa- 
tion is required to account for the very large (negative) binding energy 
of the resulting compact object (of the order of —lO^^erg) 

Indeed, it appears that the only plausible mechanism to carry away the 
bulk of the binding energy of the collapsing star, leading to a neutron 
star or black hole is neutrino emission [4]. 

• In the diffusion approximation, it is assumed that the energy flux of 
radiation (as that of thermal conduction) is proportional to the gradient 

of temperature. This assumption is in general very sensible, since the 
mean free path of particles responsible for the propagation of energy in 
stellar interiors is in general very small as compared with the typical 
length of the object. Thus, for a main sequence star as the sun, the 
mean free path of photons at the centre, is of the order of 2 cm. Also, 
the mean free path of trapped neutrinos in compact cores of densities 
about 10^^ g.cm.~^ becomes smaller than the size of the stellar core 
[5, 6]. 

Furthermore, the observational data collected from supernovae 1987A 
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indicates that the regime of radiation transport prevaihng during the 
emission process, is closer to the diffusion approximation than to the 
streaming out hmit [7]. 

Accordingly we shall restrict here to this later case, being aware that 
there are situations in stellar evolution where that approximation fails. 

During their evolution, self-gravitating objects may pass through phases 
of intense dynamical activity, with time scales of the order of magnitude of 
(or even smaller than) the hydrostatic time scale, and for which the quasi- 
static approximation is clearly not reliable, e.g., the collapse of very massive 
stars [8], and the quick collapse phase preceding neutron star formation, see 
for example [9] and references therein. In these cases it is mandatory to take 
into account terms which describe departure from equilibrium, i.e. a full 
dynamic description has to be used. 

Here we are mainly concerned with the quick collapse phase, which im- 
plies that we have to appeal to a hyperbolic theory of dissipation. The use 
of a hyperbolic theory of dissipation is further justified by the necessity of 
overcoming the difficulties inherent to parabolic theories (see references [1], 
[10]-[23] and references therein). Doing so we shall be able to give a descrip- 
tion of processes occuring before thermal relaxation. 

Many analytical solutions of Einstein's field equations with dissipative 
fiuids carrying heat fiow have been studied (see [24] for references up to 1989 
and [25], [26] for more recent ones). 

In this vein here we present some models of conformally flat dissipative 
spherical collapse with shear-free motion. We match our models to a radiat- 
ing null field described by the outgoing Vaidya spacctime. 

Although the shear free and the conformally flat conditions are introduced 
here in order to simplify calculations, it is worth noticing that these condi- 
tions generalize physical assumptions widely used in astrophysics. Indeed, 
the shear free condition in the Newtonian regime describes the homologous 
evolution and the role of shear has been extensively considered in general 
relativity [27]. On the other hand it is well known that the conformally 
flat condition implies in the perfect fluid case the homogeneity of the energy 
density distribution [2]. 

The paper is organized as follows. In section 2 the field equations, con- 
ventions and junction conditions are presented; in section 3 we present the 
general solution corresponding to the conformal flatness condition; in sec- 
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tions 4 and 5, particular analytical solutions arc given, and a specific model 
is constructed from one of them; finally a brief conclusion is presented. 



2 The fluid distribution and the interior space- 
time 

We assume a sphere of collapsing perfect fluid with heat flow. Its spher- 
ical surface S has center and is flUed with radially moving perfect fluid 
conducting heat flow, so having energy momentum tensor 

Tap ^ ilJ'+P) WaWp + pQajB + qocWj3 + t^a?/?, (l) 

where ^ and p are the proper density and pressure of the fluid, its unit 
four-velocity, the heat conduction satisfying goto" = and Qaji is the 
metric tensor of spacetime. 

We choose comoving coordinates within S and impose shear-free fluid 
motion which allows the metric be written in the form (see [28] for details) 



ds' = -A^dr + [dr' + [dO^ + sin^ 6'#^jJ , (2) 

where A and B are only functions of r and t. We number the coordinates 
— t, — r, x'^ — 9 and x^ — (f) and then we have the four-velocity given 

by 

wa = -ASl, (3) 

and the heat flows radially, 

g° = qSl (4) 
where g is a function of r and t. In these coordinates the equation of the 
boundary surface S is given by r = = constant. 

The spacetime described by (2) has the following non-null components of 
the Weyl tensor Cap-ys, 



^4 



C2323 = yS'sin^^ 




(5) 



and 



' R \ ^ / 73 \ 2 



C2323 = -r" (^^ j sin^ ^Coioi = 2t' j sin^ ^Co202 

= 2r2 (^1)' Co303 = -2r2 sin^ eC^2i2 = -2r2Ci3i3, (6) 
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where the primes stand for differentiation with respect to r and dot stands 
for differentiation with respect to t. 



2.1 The field equations 



The non null components of Einstein's field equations Ga/3 = SirTais, where 
Gai3 is the Einstein tensor and T^/^ is given by (1), with metric (2) are 

2 



Goo = --^ 



^B" B'\ AB' 
~B' " Vb] ^ r 5" 



(7) 
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B^ 
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B (B\ AB 
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sin^^ 



,B^ 
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-2 



'A 

AB 



A = -SnqAB^ 



(9) 



(10) 



The mass function m(r, t) of Cahill and McVittie [29] is obtained from 
the Riemann tensor component .^23^'^ and it is for metric (2) 

2 



m{r, t) 



(rB) 



p 23 
-rCo 



r^B 



1-23 



B' 



'B^ 



r'^B' 



(11) 



For studying the dynamical properties of the field equations and following 
Misner and Sharp [30] , let us introduce the velocity U of the collapsing fluid 
as 

U = rDtB{<0 in the case of collapse), (12) 
where the proper time derivative D^, is given by 

1 d 



^^-Adi- 



(13) 
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2.2 Junction conditions 



If the collapsing fluid lies within a spherical surface E it must be matched to a 
suitable exterior. Since heat will be leaving the fluid across S, the exterior is 
not vacuum, but the outgoing Vaidya spacetime which models the radiation 
and has metric 



2m{v) 
P 



dv^ - 2dvdp + p^{de^ + sin^ 9d(t)^), (14) 



where m.{v) is the total mass inside S and is a function of the retarded time 
V. In (14) p is a radial coordinate given in a non-comoving frame. 

The conditions for the matching of these two spacetimes (2) and (14), are 
the Darmois conditions [31], which using the field equations (8-10) and the 
mass function (11) imply [24] 

Pe = {qBh, (15) 

(rBh = PE, (17) 

(2^-Y^-^^l-™^^)' 

A.dt^^l-'-f + 2'£y\v, (19) 

where L is defined as the total luminosity of the collapsing sphere as measured 
on its surface and is given by 

and where 

dm , 
= ^ (21) 

is the total luminosity measured by an observer at rest at infinity. 
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3 Conformally flat solution 



Here we impose conformal flatness to the spacetime given by (2), i.e. all its 
Weyl tensor components must be zero valued. Then it can be shown that 
metric functions A and B take the form (see [28] for details) 



Ci(t)r2 + 1 

where Ci is an arbitrary function of t and 

1 



(22) 



B = 



C2{ty + Cs{ty 



(23) 



where C2 and C3 are arbitrary functions of t. 

Substituting solution (22) and (23) into (8), (9) and (11) we obtain, 



Snp — 



(Cir2 + 1)^ 



2{C2r^ + C^){C2r^ + C3) - 3{C2r^ + C3)' 



(24) 



-2- 



-{C2r' + Cs){C2r' + Csy 



Cir^ + 1 



Cir2 + 1 

C2(C2 - 2CiC^y + CsiCCs - 2C2 
Snq = A{CsC, - C2) ( r. 



Finally, from (15), (25) and (26) we have 

S{CW + Cs)^ Cy{C2r^ + Cs) 



C2r' + Ca - 



+2 



2 C2r2 + C73 

{cy + 1 



C2r2 + C3 



Cir2 + 1 
2 



Cir^ + 1 



- 2{CsCi - C2)r 



(25) 
(26) 



C2(C2 - 2C,Csy + CsiCCs - 2C2) 



0. (27) 



Please note that a misprint in (27) appearing in [28] has been corrected here. 

In the following sections we shall obtain some analytical solutions satis- 
fying (15)-(26). 
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4 Solution I 



In order to integrate (27) let us assume 

C2 = aCs 
Ci = Const 
a — constant 

then replacing (28) into (27) we get: 



(28) 



rr 2(Ci-a)r. • 2(Cir| + l) 



arv + 1 



2 [a(Qi-2(7i)r|+((7i-2a)]C| - 



s T -L (arl + 1) 

which, in terms of the new variable C^it) — u~'^{t), becomes: 
2(Ci - «)rs . (Cir| + 1) 



u 



ary + 1 



-u 



s -r X ^ar| + 1)^ 

The equation above allows the following three solutions: 



(29) 



[a{a - 2Ci)r|; + (Ci - 2a)]u = (30) 



Case I 



(Ci - afrl + (Cir| + l)[a{a - 2Ci)r| + {d - 2a)] > 



Cs{t) = [/3ie 
+/32e 



(Ci-a)T-s + -y/(Ci-a)2r2+(Cir2+l)[a(a_2Ci)r|+(Ci-2a)l 

( ^* 

(Ci-a)rs-^(Ci-c«)2r2+(Cir-2+l)[a(a-2Ci)r-2+(Ci-2a)] 



(31) 



Case II (Ci - a)V| + (Cir| + l)[a(a - 2Ci)r| + (Ci - 2a)] < 

Cs{t) = [e'^\p^ ^^^^ V( c.-.)2.2 +(c..2 ^i)K.-2C0.2 ^(c.-2a)i ' 



„-^^ V(gi-«)^^|; + (gi4+l)[«(«-2Ci)r|+(Ci-2a)] ^^^j_2 



(32) 



Case III (Ci - a)V| + (Cir| + l)[a{a - 2Ci)r| + (Ci - 2a)] = 

C3it) = (A + (33) 
This solution reduces to the one found in [26] when a = 0. 
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4.1 Calculation of the temperature 

It is worth calculating the temperature distribution, T{r,t), for our model, 
through the Maxwell-Cattaneo heat transport equation. For simplicity we 
shall consider here the so called "truncated" version, for which it reads [18], 
[21], 

rh'^^w'^qf,.^ + = -Kh'^^{Tp + Tap), (34) 

where r is the relaxation time, n the thermal conductivity and h""^ — g'^^ + 
w°'w^ the projector orthogonal to to". Considering (2-4) then (34) becomes 

T{qB)B + qAB"^ = -K{TAy. (35) 

In our case the integration of (35) gives 

T{t,r) = [f{t) + ^^0:^,{C,C, + Ci) 



(36) 



ar2 + l l\ Cir'2+1 



The integration function f{t) may be easily related to the central tem- 
perature Tc{t) 

m=T,{t)-^{CsCs + Ci) (37) 
then we may write for the temperature 



T{t,r)^[T,{t) + 



r (a-Ci)r-2 



4^K^^( ar^2+/)] 



ar^+1 



(38) 



where C3 is given by either (31), (32) or (33). 

In case Ci — a our system becomes a collapsing (non-dissipative) Pried- 
mann dust sphere, as it can be checked from (24)-(26). In this latter case 
the temperature, as expected, is homogeneous {T{t) = Tc{t)). Thus, mod- 
els (31)-(33) provide examples where inhomogeneity is directly related to 
dissipation. 
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5 Solution II 

Another solution, with an interesting physical interpretation may be found 
by introducing the variable v, defined as 

C2(t)r| + C3(t) = v-\t) (39) 

into equation (27), obtaining: 



(40) 



Then, defining 

Z{t) = Cir| + 1 (41) 

equation (40) becomes: 



2Zv" 



^ipl{v-^ - 2v-^Cs{l + Z) + 3ZCI) = 



(42) 



Next, introducing the new variable C{t) through C3{t) = v '^C{t), (42) 
may be written as: 



V _ Z_v I 2_v _ 2Z_V(j _|_ Z_Q 

u Z V rj2 V rj^ V »"e 

-4-(l - 2C(1 + Z) + 3ZC^) = 



(43) 



We shall integrate this last equation by assuming: 

(44) 



m = 1 

C{t) = C= const 



Observe that with this specific choice, this solution becomes a particular case 

i-c 



of solution I with Ci = Constant = and a — ^ 



Then (43) becomes 

2 . 1 



v + —(l-C)v-^(l-AC + SC^)v^O (45) 
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whose general solution is: 

v{t) ^ {ky^' + kle-^'')e^'' 

with 



c-i 

P2 = IV2-6C + 4C2 
ki and k2 are integration constants 

Then for the velocity as defined by (12) we have: 

[-2(1 - Cr- - C]r^ + [C - '-f] 



U^-r- 



(1-C)J + C 
which, evaluated at the boundary surface gives 

^ kl + A:ie-2/32* 

And for the functions C2 y C3 we obtain: 



C^{t)^v{tyC 



(fcfe'32t+fc|e-''2t)2 



'^2V''J — r|(jfc2ef2t+fcie-/32«)2 



(46) 



(47) 



(48) 



(49) 



(50) 



For this solution, the physical variables become, using (24-26) and (50): 



Hrl = 



C{l-C)-rri-[{l-C){-y-rCY 



(51) 



l-C 



((1 - C){-f - 2C) + ((1 - C){^f + C)2(^ 



r 



2 (l-C) i; 
= ^ 1 



[{l-C){-f + Cf 



+ 3C-1) 
(52) 

(53) 



which, evaluated at the boundary surface take the form: 



11 



2 .-,2 



^2 



+ C(l - C) 



(54) 
(55) 
(56) 

The expression for the temperature for this model may be obtained easily 
by putting Ci = Constant = and a = into (38), then one obtains 



2 (1-C) r^v 
Ps^s = i 

[I - C) r^v 



,6 



Tit,r) = [m + ^-^{CsCs + Cl 



+ ^ln{ar^ + l)]{ar^ + l) 



(57) 



where C3(t) is given by (50). In this case, the system becomes a collapsing 
(non-dissipative) Friedmann dust sphere, when a = (C = 1). In this latter 
case the temperature, as expected, is homogeneous {T{t) = Tc{t)). 

The last term on the right hand side of expression (57) exhibits the influ- 
ence of dissipation on the temperature, with respect to the non-dissipative 
case, as calculated from the non-causal (Landau-Eckart) [32, 33] transport 
equation, whereas the second term describes the contribution of relaxational 
effects. The relevance of such effects have been brought out in recent works 
(see [25] and references therein). In particular it is worth noticing the increas- 
ing of the spatial inhomogeneity of temperature produced by the relaxational 
term, an effect which has been established before [34]. 

Let us now present a very simple model based on the solution above. The 
purpose here is not the modelling of any specifical astrophysical scenario, 
but rather to show the feasibility of these solutions as starting point for such 
modelling. 

Thus, let us consider the following choice of constants and initial values: 



and 



k-^ ~\~ — 1 
C = 1 + 10"^ 

?7s(0) = -2.5 X 10" 



(58) 
(59) 

(60) 
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Then it follows from (47) and (49) that: 

A - — (61) 
ft = 1-4142 . 10-3(1 + 10-^) 

kj = .0577 (63) 

kl = .9423 (64) 

With these values, it follows at once from (54)-(56) that ^s(O) > qsiO) — 
Ps(0) > 0. Furthermore since /x is a decreasing function of r, these inequal- 
ities hold for all points within the sphere. As time goes on, the velocity of 
the boundary surface decreases (in absolute value), eventually changing of 
sign for a finite time value. However physical variables remains acceptable 
(in this model) for values of the boundary velocity close to zero, but still 
negative. 

Thus our model describes an initially contracting and radiating sphere, 
approaching the equilibrium. For later times there is a bouncing of the 
boundary surface, however, physical variables become unphysical for this 
values of t, and the model is restricted to the collapsing regime only. 



6 Conclusion 

We have presented some exact analytical solutions to the Einstein equations, 
describing spherical dissipative shear-free and conformally flat collapse. The 
solutions are matched to the outgoing Vaidya radiating spacetime. Besides 
their simplicity, the merit of the models resides in the fact that they exhibit 
in a very clear way the influence of relaxational effects on the temperature, 
and thereby on the evolution of the system. 

In this respect we would like to stress the modifications in the temperature 
profile of the models, produced by the relaxational efects. This fact cannot 
be over emphasized. Indeed, different temperature profiles, are not only 
associated with different patterns of evolution, but also, affect the luminosity 
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profile, which is the most important element of observational evidence in the 
study of dissipative collapse. 

In the same line of arguments, it is worth noticing that the resulting 
temperature profile for each model, will depend on the specific theory of 
transport employed in its calculation. Therefore, such models might be used 
as test-bed for different rclativistic theories of dissipation. 

The specific example presented at the end of the previous section, mod- 
els a dissipative collapsing configuration approaching equilibrium, with all 
physical variables exhibiting appropriate behaviour. This support further 
our believe that the presented solutions may be suitable for describing astro- 
physical scenarios involving dissipative collapsing objects. 

Finally, it is also worth noticing that density inhomogeneities are directly 
related to dissipation, while the space-time remains conformally flat. In the 
non-dissipative limit {q = L ~ 0), all models become homogeneous dust 
balls matched to Schwarzschild spacetime. 

This reinforces doubts ( see for example [35] and references therein) on the 
proposal that the Weyl tensor [36] or some functions of it [37], could provide 
a gravitational arrow of time. The rationale behind this idea being that 
tidal forces tend to make the gravitating fluid more inhomogeneous as the 
evolution proceeds, thereby indicating the sense of time. However, as shown 
in [2], density inhomogeneity, besides Weyl tensor (and the anisotropy of 
pressure), also depends on dissipation. The solutions obtained here, clearly 
bring out that dependence. 
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